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While the hehcal character of the edge channels responsible for charge transport in the 
quantum spin Hall regime of a two-dimensional topological insulator is by now well estab- 
lished, an experimental confirmation that the transport in the edge channels is spin-polarized 
is still outstanding. We report experiments on nanostructures fabricated from HgTe quantum 
wells with an inverted band structure, in which a split gate technique allows us to combine 
both quantum spin Hall and metallic spin Hall transport in a single device. In these de- 
vices, the quantum spin Hall effect can be used as a spin current injector and detector for 
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Q the metallic spin Hall effect, and vice versa, allowing for an all-electrical detection of spin 

r^ polarization. 

O 

> I. INTRODUCTION 

X 

C3 The discovery that HgTe quantum wells (QWs) with an inverted band structure are 2- 

dimensional topological insulators has generated a great interest in this novel state of quantum 
mattei^ ^. When the thickness d of the HgTe QW is increased beyond a critical value dc, a quan- 
tum phase transition turns a conventional insulator into its topologically non-trivial counterpart. In 
this so-called quantum spin Hall (QSH) phase"^^, current-carrying states are confined at the edge 
of the sample, while the bulk is insulating. These edge states are protected against backscatter- 



ing from non-magnetic impurities^ ^ and their propagation direction is helical, i.e. that opposite 
spin states counter-propagate along a given edge of the sample. When the applied gate voltage 
places the Fermi level inside the bulk gap, two-terminal transport experiments measure a quan- 
tized conductance of 2e^//i with e the electron charge and h the Planck constant, independent 
of the sample width, which constitutes strong evidence for edge state conduction^. More recent 
nonlocal transport measurements in the QSH regime unambiguously establish that transport oc- 
curs through topologically protected edge channel^^^. While edge state conduction in the QSH 
regime is thus experimentally well established, there exists so far no direct experimental evidence 
that the transport in the helical edge states of 2D topological insulators is spin-polarized, which is 
a fundamental characteristic of this new state of matter. 

In this work, we construct novel devices (Figs. 1 and 2a) that enable us to study the spin 
polarization of the QSH edge states by purely electrical means. First of all, these devices allow us 
to detect the spin polarization of the QSH edge states (Fig. lb) via the inverse spin Hall effecF^^ 
(SHE^). Second, our devices enable us to show that because of their helical nature, the QSH 
edge states can be used as a detector of spin current (Fig. la). In our devices, the spin current 
is generated by the intrinsic ballistic spin Hall effecP^ (SHE) exhibited by a HgTe QW in the 
metallic regimeP'. These two experiments establish for the first time spin polarization of the helical 
edge states in topological insulators, and also demonstrate potential applications of the QSH effect 
for spintronic devices. 

II, PRINCIPLE OF THE EXPERIMENT 

Before presenting our results, we first describe the principle of our experiment in more detail. 
Since the magnetic field originating from spin polarized carriers in helical edge channels is too 
small to be detected directly, we have designed an experiment that converts magnetic information 
into an electrical signal. Figure 1 illustrates the idea of the experiments, which are performed on 
an H-shaped mesa structure (which we call 'H-bar') in which the carrier concentration in the two 
legs of the 'H' can be adjusted separately. Consider the situation illustrated in Fig. 1, where the 
bottom leg is metallic (indicated by the green color, and either n- or p-type) and the top leg is 
tuned into the QSH regime (indicated by the yellow color), with the counter-propagating helical 
edge channels depicted as blue and red trajectories. We perform two separate complementary 
experiments. 



In Fig. la, the current is injected into the metallic part of the structure (contacts 3 and 4) 
while a voltage signal is detected across the top leg (contacts 1 and 2), which is gated into the 
QSH insulator state. The inverted band structure in HgTe results in a large spin-orbit couplin^^^^, 
which has previously enabled us to observe a ballistic intrinsic SHE in a small H-bar structure with 
a homogeneous carrier profile^^. Similarly, when in the experiment of Fig. la a charge current is 
injected into the metallic leg, the intrinsic SHE will induce a separation of carriers with opposite 
spin polarizations toward opposite edges of this leg. This leads to a difference in chemical potential 
for opposite spin states in the area where the metallic part of the structures borders on the QSH 
region. The spin polarized helical edge channels coming from the QSH region couple selectively 
to the chemical potential of the matching spin species in the metallic region and transfer this 
difference in potential to voltage contacts 1 and 2. For non-spin-selective edge channels the voltage 
signal is expected to be zero, while for the spin polarized QSH edge channels we expect a nonzero 
signal. Thus the observation of a nonlocal signal in this configuration is evidence that the metallic 
leg develops an intrinsic SHE, as well as that the helical edge channels are spin polarized in the 
QSH insulator regime. 

In the reverse configuration of Fig. lb, the current is injected (contacts 1 and 2) into the area 
of the sample that is gated into the QSH regime, while a nonlocal voltage drop is measured across 
the metallic leg (contacts 3 and 4). In this configuration, the spin polarized helical edge channels 
inject a spin polarized current into the metallic leg, causing a local imbalance in the chemical 
potential of spin-up and spin-down polarized carriers. Due to the SHE^ (see Refs.'^^^^, the spin 
current in the metallic region induces a voltage between contacts 3 and 4. Again, this voltage can 
only develop provided the helical edge channels are spin polarized, and the metalhc leg exhibits 
the SHE 1. 

A possible complication in both of the above experiments is the detection of a stray spreading 
voltage. In the configuration of Fig. la, this could result from a voltage drop in the metallic leg 
along the interface with the area in the QSH insulator regime, while in Fig. lb, the finite distance 
between in- and outgoing edge channels at this interface could produce a similar effect. However, 
in practice such stray voltages are strongly reduced by the exact layout of the experiment, the 
quasi-ballistic nature of the transport in the metallic leg and the finite width of the edge channels. 
In the supplementary online material (Section VI), we present theoretical models and experimental 
evidence showing that any residual stray voltages can be neglected as compared with the actual 
signals obtained in the experiments described below. 



Ill, SAMPLE FABRICATION AND TRANSPORT CHARACTERIZATION MEASUREMENTS 

Our H-bar structures are fabricated from inverted HgTe/HgCdTe type-III QWs with a nom- 
inal well width of 9 nm, located 74 nm below the surface. At a temperature of 4 K (and for 
a grounded gate), the carrier density is n ^ 4 x 10^^ cm^^ . The carrier mobility is then 
/i ^ 1.1 X 10^ cm^/(V-s), yielding an elastic mean free path larger than 2 /im. The devices 
are patterned using optical and electron beam lithography, with dimensions as indicated in Fig. 1. 
In order to control the carrier density, the device is gated by Au gate electrodes separated from the 
sample surface by a 110 nm thick insulating Si3N4/Si02 multi-layer stack. By applying a voltage 
Vgate to the top gates, the electron carrier density of the QW can be adjusted, going from an n-type 
behavior for Vgatc > through the bulk insulator state into a p-type regime for Vgatc < 0. For 
reasons of comparison, the experimental data in Fig. 2, 3 and 4 are plotted as a function of a nor- 
malized gate voltage V*^^^ = legate — Khr, where the threshold voltage Vj^hr is defined as the voltage 
for which the resistance is largest in a particular fixed reference measurement. As is evident from 
the characterization data in Fig. 2b and c, which were obtained from a Hall bar fabricated from the 
same wafer material as the H-bar nano structures, we find that for gate voltages V*^.^ > 0.5 V the 
QW is n-type metallic, and for legate ^ —0.5 V it is p-type metallic. The split-gate design (gates 
1 and 2) of Fig. 1 provides an independent control of the carrier density for each leg of the H-bar 
structure, enabling us to gate one part of the sample into the QSH insulator regime and the other 
part into either n- or p-type metallic regimes. An electron microscope picture of the actual device 
structure is shown in Fig. 2a. The transport measurements are done at a constant temperature of 
1.8 K employing quasi-dc low frequency (13 Hz) lock- in techniques using a voltage bias below 
100 yuV. 

IV. EXPERIMENTAL RESULTS 

While experiments have been performed on a variety of different devices and yield similar re- 
sults, for reasons of consistency we will discuss here a single device with dimensions as indicated 
in Fig. 1. The results of the experiments are shown in Figs. 3 and 4, corresponding to the mea- 
surement configurations of Figs, la and b, respectively. In the upper three panels of the figures, 
the nonlocal resistance is plotted as a function of gate 1, while in the lower panel, gate 2 is swept. 

Figure 3 corresponds to the layout of Fig. la, and the detected nonlocal signal can consequently 



be denoted as -R34,i2, i-e- the voltage measured between contacts 1 and 2 divided by the current 
passed between contacts 3 and 4. When we sweep the gate on the injector area (gate 2) while 
the detector is tuned into the QSH regime (V*^^^^ = 0), we observe (lower panel in Fig. 3) a pro- 
nounced maximum around V*^^^^ = 0, and smaller but finite values on both sides. The signal 
around V^g^tci = reaches approximately the quantized value observed in our previous experi- 
ments on nonlocal transport in the QSH regimeP. We attribute the slight deviation from perfect 
quantization to imperfect gating in the not gate-covered region between gates 1 and 2. Imperfectly 
gated regions in the sample can act as dephasing centers for edge electrons, which can lead to a 
deviation from the expected quantized nonlocal resistanceP^. In addition, in HgTe QW devices 
subsequent gate voltage sweeps can charge interface trap states in a different wajP^, leading to 
different dephasing effects and a different magnitude of the deviation from quantized resistance 
for each gate voltage sweep. 

Apart from the large signal in the QSH regime, the measurements also exhibit a non- vanishing 
nonlocal signal when the area underneath gate 2 is metallic, either n- or p-type, and thus corre- 
sponds to the injector region depicted in Fig. la. The origin of this finite signal becomes more ev- 
ident when the injector gate voltage is set at a fixed value either in the p-type (V*^^.^2 = —0.75 V < 
0) or in the n-type metallic regime (V^*ate2 = 1.0 V > 0) while the voltage on gate 1 is swept 
(top panel of Fig. 3). Evidently, a significant increase in the nonlocal signal is observed, with a 
peak when the detector is exactly in the QSH insulator regime. This is the observation anticipated 
above: one may expect a nonlocal signal of this amplitude only when the metallic leg exhibits a 
SHE and the edge channels in the QSH leg are spin polarized. Our data also show that the nonlocal 
signal for the p-type injector (V^*ate2 = —0.75 V) is more than ten times larger than that for the 
n-type injector (V^*atc2 = 10 V). This is consistent with our experimental observations on the SHE 
signal in all-metallic HgTe QW^, where the nonlocal signal is about an order of magnitude larger 

in the p-regime than in the n-regime and results from enhanced spin-orbit splitting in the valence 
bancPni. 

Our data for the reverse configuration of Fig. lb are shown in Fig. 4. The sweep of gate 2 in 
the bottom panel now corresponds to the detection leg, and one can directly see that also in this 
configuration we observe a finite nonlocal signal (in this case -Ri2,34), even when the detector is 
metallic (red and green arrows). The upper panel shows the effect of sweeping the injector leg (gate 
1), and indicates that the nonlocal signal peaks when the injector is in the QSH state. As in the 
previous configuration of Fig. 3, we observe an order of magnitude increase in the nonlocal signal 



when the metallic detector is p-type (V^*ato2 = —0.82 V) as compared with an n-type detector 
(K*ate2 = 1-2 V). As noted above, our observation of the nonlocal signal is evidence that the 
helical edge channels generate a spin accumulation at the interface between the QSH injector and 
the metallic detector, which responds by the SHE^. 

The results in Figs. 3 and 4 look very similar and, in fact, are expected to do so on account 
of the Onsager-Casimir symmetry relations for the nonlocal resistances Rmn,ki in a four-probe 
devicdi^, 

Rmn,kl(^) ~ Rkl,mn{ — ^), (1) 

where the first pair of indices refers to the current probes, the second pair refers to the voltage 
probes, and B is the magnetic field. In the present setup, the magnetic field is zero and we expect 
-^34,12 = -Ri2,34- One possible explanation for the small deviations from exact Onsager-Casimir 
symmetry observed in Fig. 3 and 4 is the random charging effects of pinned inhomogeneities (or 
'trap states') mentioned earlier. Two subsequent gate voltage sweeps can result in a different in- 
terface potential due to these charging effects^, which changes the internal state of the conductor. 
Note however that the symmetry between Fig. 3 and 4 is more accurate in the doped regimes away 
from the nominally insulating regime, which is expected since a higher carrier density can more 
effectively screen the interface trap potentials and thus make the internal state of the conductor 
less sensitive to trap charging effects. 

V. MODELING AND DISCUSSION 

In order to better understand the experimental results, we have performed semi-classical Monte 
Carlo calculations to obtain a theoretical estimate of the nonlocal resistance based on the sample 
geometry (Fig. 1). We focus on the setup illustrated in Fig. lb, where the QSH insulator acts as 
a spin injector and the metallic region detects the spin polarization of edge channels through the 
SHE^^ We calculate the nonlocal resistance -Ri2,34 when the current is driven between contacts 
1 and 2 and the voltage is measured between contacts 3 and 4. -Ri2,34 can be expressed in terms 
of the transmission coefficient^^ Tij for the metallic region only (Eq. SI of the supplementary 
online material). The Tij are calculated within the semiclassical Monte Carlo method^", which is 
a reasonable approximation for Fermi wavelengths Xp <^ L where L ~ 1 fim is the characteris- 
tic linear size of the device (Fig. 1 and 2a). Electrons are injected at the QSH-SHE^ interface 



(yellow-green interface in Fig. lb), and propagate quasi-ballistically into the metallic T-structure 
(green region in Fig. lb) according to semiclassical equations of motion^'. These equations are 
derived using an effective 4-band model for HgTe Qw^ which explicitly contains the effects of 
intrinsic spin-orbit coupling due to atomic coupling between bandsSEl x^is intrinsic spin-orbit 
coupling can be visualized as resulting from a Rashba field due to the edges of a typical mesa 
structure used in experiments'"*. Details of the calculation are included in the supplementary on- 
line material. 

We find that the conversion of the spin signal to the electrical signal through the SHE ' is 
dominated by the intrinsic spin-orbit interaction while stray contributions due to voltage spreading 
are negligible (see Fig. S3 in the supplementary online material). Fig. 5 shows the theoretical 
prediction of the nonlocal resistance signal as a function of the carrier concentration in the metallic 
detector. (Note that the semiclassical simulation breaks down when the chemical potential is too 
close to the insulating gap.) The scattering induced by the intrinsic spin-orbit interaction is more 
effective when carriers have smaller kinetic energy, and therefore smaller wave vectors at the 
Fermi leveP^. Since the effective mass in the p-regime is larger than that in the n-regime^^, for 
comparable densities the kinetic energy will be smaller in the p-regime. This can explain the 
larger nonlocal resistance signal for the p-regime in comparison with the n-regime, as well as the 
decrease of the signal upon increase in carrier concentration. 

VI. CONCLUSION 

In conclusion, we have presented nonlocal transport experiments on split-gate HgTe QW which 
establish for the first time the spin polarization of the helical edge states in topological insulators. 
Together with our previous experiments on conductance quantizatioiP and nonlocal edge state 
transporr', the present data constitute the final piece of evidence needed to fully confirm the orig- 
inal predictions for this novel state of mattei' ^^. Our experiments yield a robust signal, in both 
experimental configurations where the helical edge channels are used for spin injection and detec- 
tion. This opens a novel route for spintronics experiments in two-dimensional electron systems at 
zero external field. 
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VIII. SUPPLEMENTARY INFORMATION FOR ARTICLE ' SPIN POLARIZATION OF THE 

QUANTUM SPIN HALL EDGE STATES' 

In this supplementary online material, we present further theoretical results and details of our 
semiclassical Monte Carlo calculations, as well as additional experimental data for the nonlocal 
resistance in structures with different injector size. These establish that the helical property of the 
quantum spin Hall edge states can be detected via the inverse spin Hall effect, and that the spin 
Hall effect can be converted into a charge signal through the helical edge channels of the quantum 
spin Hall state. 

IX. INTRODUCTION 

This document provides the details of our theoretical approach for the description of the hybrid 
quantum spin Hall - metallic spin Hall (QSH-SHE) 'H-bar' structures illustrated in Fig. 1. of 
the main text. More specifically, we calculate the nonlocal four-terminal resistance -Ri2,34 in the 



configuration illustrated in Fig. lb and Fig. S6 where the quantum spin Hall (QSH) insulator 



state acts as a spin current injector, and the metallic region as a spin current detector through the 
inverse spin Hall effect (SHE^^). The nonlocal resistance -R34,i2 in the reverse geometry (Fig. la), 
with voltage and current probes exchanged, is the same as -Ri2,34 by virtue of the Onsager-Casimir 



symmetry relations in the absence of a magnetic field (Refs. [181 [I9|and Eq. (1) of the main text). 
We find a good qualitative agreement between the experimental and theoretical values (Fig. 4 and 
5 of the main text). 

The main idea of our approach is the following. We consider the H-bar geometry illustrated in 



Fig. S6 The yellow region on the left is in the QSH regime, while the green region on the right is 



in the metallic regime, either n-type or p-type. Figure S6 represents the same setup as in Fig. lb, 
where LI corresponds to the length of the part of the structure in the metallic regime. 

We calculate the four-terminal nonlocal resistance -Ri2,34 = (V3 — V4)// where I = Ii = —I2 
is the current injected into lead 1 and collected in lead 2, and V3 — V4 is the voltage difference 
between leads 3 and 4. An expression for -Ri2,34 in terms of the transmission coefficients Tij from 
lead j to lead i,i,j = 1, . . . , 4, can be derived from the multiprobe Landauer-Biittiker formulcP^^, 
with Tij = Tji from time-reversal symmetry'22l(TRS). The coefficients Tu = T21 = 1 in the QSH 
region are universal and are obtained simply by counting the number of edge states leading from 
one contact to the other (see Ref. 9). In the same spirit, the coefficients T34 = T43 can be 
approximated by counting the number of channels^^ N^. ~ kpWl-K (per spin) in the vertical arm 



of width W = W2 in the metallic T- structure (Fig. S6), with kp the Fermi wave vector in the 
metallic region. The nontrivial part of the calculation is to determine the coefficients T13, T14, T23 
and T24 for the transmission between the contacts attached to the QSH region and those attached 
to the metallic region. Since the right T-structure is metallic, for high enough densities the Fermi 
wavelength Xp <^W is much smaller than the dimensions W of the T-structure and a semiclassical 
approach becomes justified. In other words, there is a large number of channels Nc ^ 1 in the 
metallic region and we can neglect the quantization of motion in the transverse direction. 

Our approach for the calculation of T13, T14, T23, T24 is based on the Monte Carlo method for 
the simulation of semiclassical electron transport in semiconductors^S. This approach has been 
used successfully for the study of magnetotransport in multiprobe conductors^SEHlI]^ ^^^ ^^ partic- 
ular for the study of quasiballistic transport in HgTe quantum well^. The procedure consists in 
injecting electrons with well-defined positions and momenta at the left of the metallic T-structure; 
the QSH edge states then act as 'injectors' of electrons into the metallic region. The electrons 
then propagate into the ballistic region according to semiclassical equations of motion^, which 
include elastic scattering on the geometric boundaries of the T-structure, elastic scattering on non- 
magnetic impurities through the inclusion of a phenomenological momentum relaxation time, and 
more importantly, the effect of spin-orbit coupling through the inclusion of a Berry phase tenrP^ 
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which acts as a magnetic field in momentum spac d^'^^^ *^. This last term is responsible for the SHE 
and SHE^. One then simply counts the fraction of electrons which reach contacts 3 and 4, from 
which the classical transmission probabilities T13, T14, T23, T24 can be extracted^-. 

This document is structured as follows. In Sec.|X} we apply the multiprobe Landauer-Biittiker 
formulcT^ to the geometry of Fig. S6 and give an explicit expression for the nonlocal resistance 



-Ri2,34 in terms of calculable quantities. In Sec. XI we apply the formalism of Ref.l21lto derive the 



form of the equations of motion describing carrier propagation in the metallic region, including 



the important Berry phase term. In Sec. XII, we give the details of the Monte Carlo algorithm. In 



Sec. XIII, we present and discuss our numerical results. Finally, in Sec. XIV we present additional 
experimental and theoretical results which provide strong evidence that the observed nonlocal 
signal is a direct consequence of spin-orbit coupling in our samples. 



X. MULTIPROBE LANDAUER-BUTTIKER FORMALISM AND NONLOCAL RESISTANCE 



For a four-probe phase-coherent device, the nonlocal resistance -Ri2,34 is one of several possible 
four-terminal resistances, whose expressions in terms of transmission coefficients were all worked 
out long ago by Biittikei^i^. He obtained 

V3 — V4 h T31T42 — T32T41 



-Rl2,34 = 

at zero temperature, where I = Ii 
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—I2 is the current injected at contact 1 and collected at 
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where Tjj = Ti^j is the transmission probability from lead j to lead i at the Fermi level, and 
Tij = Tji by TRS in the absence of a magnetic fielcP^. 



From Fig. S6 one can read off T12 = 1 + T[2, where 1 is universal and comes from the 



edge state propagating directly from contact 1 to contact 2 along the left edge of the QSH T- 
structure, while T^'2 i^ nonuniversal and is the probability of an electron travelling from contact 
1 along the top-right edge state of the QSH T-structure, entering the metallic region, propagating 
inside the metallic region, returning inside the QSH T-structure and propagating to contact 2 via 
the bottom-right edge state. In fact, T^g corresponds to the probability of interedge tunneling, 
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which is negligibly small for a wide enough device. Indeed, Ref. 'BT finds that for a device width 
Wl ~ 1 yum (we have Wl = 0.77 /im, see Fig. 1), the gap A opened in the edge state dispersion by 
interedge tunneling is negligibly small, A ~ 10^^ meV. Also, nonlocal resistance measurements 
(similar to those of Ref. 9 performed on the actual devices when they are entirely gated into the 
QSH regime yield the values expected from unperturbed nonlocal edge state transport (see Ref. 9 
and Fig. 2b), suggesting that interedge tunneling is negligible. Finally, we have performed fully 
quantum-mechanical, numerical calculations of the S-matrix of a QSH/metal interface in a strip 
geometry using the tight-binding version of the four-band effective model for the QSH state in 
HgTe quantum wellP, and confirm that T^g is negligible for the sample widths considered here. 
Therefore we take T[2 = and T12 = 1. 

The coefficient T34 for transmission from contact 3 to contact 4 through the metallic region 



(Fig. S6) is obtained from the semiclassical Monte Carlo calculation as follows. We first calculate 
T^4 and R^ defined as the fraction of electrons of spin a =t, \. injected from contact 3 that reach 
contact 4, or that are reflected back into contact 3, respectively. They satisfy R^ + T^ = 1 since 
the metal/QSH interface is modeled as a perfectly reflecting interface and the z component of 



the spin is conserved in our model (the electron spin is discussed in Sec. XI and XII). However, 
these coefficients assume only one transport channel (since they sum up to one^^), and furthermore 
neglect the probability of being transmitted into the QSH region T31 + T32. The metallic leads have 
Nc channels per spin at the Fermi level^ with N^. ~ A;^iy/7r and W = W2 the width of the lead 



(Fig. S6 1. Taking this into account as well as the probabilities T31, T32, the actual transmission and 
reflection coefficients T^^, R^ should satisf))^!] {T^^ + TQ + T^^ + Rl = N^. We can thus construct 
transmission and reflection coefficients satisfying this constraint by defining 

Rl = Rl[N,-{T^,+T^,)l (S2) 

T^, = f^,[N,-{T^, + TQ]. (S3) 

The total transmission coefficient is then given by T34 = X]cr=t 1 -^34- 

In the Monte Carlo procedure (to be detailed in Sec. |XII| ), the semiclassical equations of mo- 
tion, which are first order in time, are integrated numerically, starting from the initial positions 
and momenta of the charge carriers. In order to obtain the transmission coefficients, we need to 
perform two separate Monte Carlo calculations. A first calculation, where Fermi surface electrons 
are injected at the QSH/metal interface and collected at the contacts 3 and 4, yields T13, T14, T23 
and T24. The problem of the interface between a QSH insulator and a normal metal is nontrivial 
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and has been studied in Ref.|32] using quantum- mechanical scattering theory. In the present work, 
we use a simpler semiclassical approach. We take for initial conditions at the left of the metallic 
region a spatially uniform distribution, a fixed wave vector amplitude kp = |k| equal to the Fermi 
wave vector, and an angular distribution given h'}p^^ 



where 9 is the angle between the carrier wave vector k and the the positive x axis (Fig. S6), and 



j^^,^d9 P{9) = 1. The reason we choose a spatially uniform distribution is that although the 
electrons are injected from the QSH side with wave functions localized along the edge, as soon as 
they enter the metallic region their wave functions merge into the bulk and their localization length 
diverges^. A second calculation consists in injecting electrons from contact 3 and collecting them 



at contacts 4 or 3, and yields T34 according to Eq. (S3 1. The angular distribution is again given by 
P{9) = I cos 6', — I < 6* < |, but now 9 is the angle between the carrier wave vector k and the the 
negative y axis (Fig.|S6]). Positionwise, the carriers are again injected uniformly across the width 



W2 of contact 3 (Fig. S6 1 



XL SEMICLASSICAL EQUATIONS OF MOTION AND BERRY PHASE 

In this section we derive the equations of motion which contain the reciprocal magnetic field 
in momentum space, or Berry phase tenrP^'. We neglect the bulk inversion asymmetry terms 
which are a small perturbatioiP. We also neglect interband transitions and study separately the 
semiclassical dynamics of Kramers partners at the Fermi level in the metallic regime. For either 
of the n- or p-type regimes we have two degenerate bands that are related by TRS, which means 
that in our simulation we track the position and momentum of two species of particles denoted by 
t and I. The energy of the degenerate conduction bands (n-type) is E^ and that of the degenerate 
valence bands (p-type) is E_ . 

As before, we use the convention of Ref.lBTIfor the Hamiltonian. The spectrum consists of two 
energy eigenvalues, 

E±{P)=e{P)±d{P), (S4) 

where e{k^) = ~Dk^, d{e) = ^/A^WTMHW), and M'^{k'^) = M - Bk^. Each eigenvalue is 
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two-fold degenerate, with eigenstates 



f{k)) = f i<(k)) = «(k)); I - i) + {uf{k)); I - 1). 



The eigenspinors are given by 

±fk) = f ^''^^'^^^l ^ = 1 ( ±MK - iky, 

and are orthonormal, where g±{k'^) = d{k'^) =F M{k'^). We can now calculate the Berry curvatures 
for each bancpil, 



^J-(k) = -2Im/^"^('^) ^"^('^^ 



^apy'-)- ---\ Qj^ 



dkf- 



(S5) 



with a =t, i and a, /3 = x, y. We however immediately observe that due to TRS, we have 

therefore we only need to calculate the Berry curvature for spin f. Furthermore, Vtpa = —^ap is 
antisymmetric from the definition Eq. ( [S5| ). In two dimensions, this antisymmetric tensor has a 
single component Vt^y and we can define a pseudoscalar 

fi±(k) ^ ie"^^fij^(k) = fi^^^k), 

where e""^"^ is the Levi-Civita symbol and a = z necessarily. We obtain 



2 d(fc2) [A2F^^±(P)M(A;2)]2- 

Note that the physical units of the quantities A, fc, f2 are given by [A\ = eV-A and [A;] =A^^, hence 
[Q] =A2. The anomalous or Hall velocities are given by^ 

Mx± = -hn^fk^,^ = -am^k^,, 
hSyt = -hn^fkt^ = <yhn^k%„, 

with a = ±1. Since the bands are doubly degenerate, the normal velocity is independent of spin 
and is simply given by 

hv± = a±(k^)k with a±(k^) = -2D ± ttt^t^^, 

d[k^) 
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where [a] = eV-A^. The semiclassical equations of motion finally take the form 



hx^ = a±((k^)2)A;J, - ann±((k±)2)fcj^, (S6) 

hy^ = «±((k^)2)fcj, + an(7±((k±)^)fcj,, (S7) 

ni±, = F-i'(r^,k^,t), (S8) 

nfc,±, = F-i'(r±,k±,t), (S9) 

where [hr] = eV-A, and F'^°^^ is the force exerted on the particles due to collisions with the 
geometric boundaries of the sample and with impurities in the sample. The specific form of this 
term is detailed in Sec. IXIII 

XII. SEMICLASSICAL MONTE CARLO ALGORITHM 

For simplicity, we assume that the probability of reflecting into a state with opposite spin is 
very small, so that we are effectively simulating the semiclassical dynamics of a two-component 
gas, with the two components evolving in a perfectly independent manner. 

The collision force is given by F'^"" = F'^'^ + F'™p where F^*^ is the force due to collision on the 
sample boundary (denoted by OS), and F'™p is the force due to collisions on impurities. The effect 
of F^'^ on a particle is implemented into the simulation as follows: by energy and momentum 
conservation, when a particle hits a boundary we simply flip the sign of the component of its 
momentum normal to the boundary (specular reflection). On the other hand, to take the effect of 
pimp jj^j-Q account we proceed as follows. We first generate a random free flight time^"^ tfree from 
an exponential distribution 

P(tfree) = -6-''-^^^ 
T 

where r is a phenomenological collision time. We expect the metallic region to be in the quasi- 
ballistic regimepland thus consider that collisions are dominated by boundary scattering. There- 
fore, we choose r > tqs where tqs ~ L/vp is the boundary scattering time with L the character- 
istic linear size of the device and vp the Fermi velocity. The impurity-free equations of motion, 



i.e. Eqs (S6)-(S9) with F'™p = 0, are then solved numerically for a time tfree- At the end of the 
free flight time, we randomize the momentum of the particle (only the direction k = k/|k| as the 
magnitude |k| = fci;' is fixed by energy conservation) according to an uniform angular distribution 
between and 2tt. This is meant to simulate isotropic scattering from rotationally invariant impu- 
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rities. Finally, a new random free flight time is generated, and the procedure starts again until all 
particles have exited the device through either contact 3 or 4. 

Since Kk^ = apart from boundary and impurity scattering, the Berry phase term vanishes 
for free propagation inside the boundaries and the semiclassical trajectories between collisions 
are simply straight lines. The Berry phase term generates a shift of position upon scattering on 
boundaries and impurities, an effect similar to the side-jump effect^^in the anomalous Hall effect, 
with the exception that here the spin-orbit coupling is intrinsic (arising from the bandstructure) 
and not the spin-orbit coupling arising from impurity potentials. A scattering event resulting in a 
change of momentum Ak produces a position shift Ar given by 

Ax^ = -afi^((k^)2)AA:J,, 
Ay^ = afi±((k^)^)AA;t, 

which can also be written 

Ar± = -af2±((k±)2)xAk±, (SIO) 



where f2^ = Q,^z and Ak^ depends on the boundary. Equation (SIO) makes explicit the interpre- 
tation of the Berry curvature f2^ as a magnetic field in momentum space. Since Ak^ is normal 
to the boundary for specular reflection, the position shift Ar^ will be along the tangent to the 
boundary. For free propagation inside the boundaries, we have 

Ar^ = ai((k^f)k^^, 

where [^] = eV^^, and At is the time between two collision events (impurity or boundary). 
Furthermore, since we are doing a Monte Carlo simulation we need to average over a large number 
of particles. We find that we get reasonably good statistics (error bars not too large) for ~ 10^ 
particles. 

XIII. NUMERICAL RESULTS 

We have performed the calculations for HgTe/HgCdTe quantum wellP of thickness d = 



89.9 A and device size (see Fig. |S6|) LI = 0.2 /im, Wl = 0.77 /im, W2 = 0.8 /xm and 
L2 = 1.665 A (such that the distance between contacts 3 and 4 is 4.1 /im). These dimensions 
correspond to those of the device discussed in the main text (Fig. 1). Note that although the 
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length of the middle segment of the H-bar is 0.7 /xm (Fig. 1), the actual length LI of the QSH 
injector region (yellow region in Fig. 1) is estimated as 0.2 fim. To avoid unphysical geometric 
resonances^SISl]^ ^g consider that the comers of the T-structure are rounded, with radius of curva- 



ture R = 100 A. For the devices studied in this work, the Dirac mass M (see Sec. XI) is estimated 
as |M| = — 6 meV due to the small gap between the HI and H2 subbands which are the lowest 
energy subbands for the quantum well thicknesses considered. 

The results are plotted in Fig. 5 of the main text. The calculated nonlocal signal is in good 
qualitative agreement with the experimental results of Fig. 3 and 4, with -Ri2,34 ~ 10^ fi in the 
p-type regime and -Ri2,34 ~ 10 fi in the n-type regime. Since we are performing a semiclassical 
simulation for the detector region, we cannot simulate the transition through the insulating gap, 
since as the Fermi level approaches the gap, the density reaches a point where our semiclassical 
approximation kpyV/n ^ 1 breaks down. Although the qualitative agreement with experiment is 
good, we observe that the experimental signal is larger. We expect that the discrepancy is due to 
the additional contribution from Rashba spin-orbit coupling^, which is not taken into account in 
the simple semiclassical approach with Sz conservation used in this work. 

The increase of the nonlocal resistance signal with decreasing density (Fig. 5) is evidence that 



the signal is generated by the spin-orbit interaction, i.e. the Berry curvature term (Sec. XI). This 



can be understood if one remembers that the Berry curvature can be treated approximately as the 



magnetic field in reciprocal space of a magnetic monopolePESl centered at k = 0. As the density is 
reduced (decreasing k), the particles are closer to the monopole and feel a stronger Berry magnetic 
field. Furthermore, due to the larger effective mass for holes (p-type) than for electrons (n-type) 
in our structures^, for comparable densities the Fermi wave vector will be smaller for holes than 
for electrons, yielding a larger Berry phase effect in the p-regime as compared to the n-regime. 

XIV. ADDITIONAL EXPERIMENTAL AND THEORETICAL EVIDENCE FOR SPIN-ORBIT 

ORIGIN OF OBSERVED EFFECTS 

We now present additional experimental and theoretical evidence that the observed nonlocal 



resistance signal is due to the spin-orbit interaction and not by some spurious effects. In Fig. S7 



we plot the nonlocal resistance as a function of gate voltage for shorter (Fig. S7 1, 200 nm) and 



longer (Fig. S7 ?, 400 nm) injector sizes. The sample parameters are the same as in the main 



text. Although the data in Fig. S7 corresponds to a metallic (p-type) injector, for V*^^ ~ in the 
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detector (QSH insulating regime) the configuration is equivalent to that illustrated in Fig. |S6| by 
Onsager-Casimir reciprocity (see Eq. (1) in the main text). In our calculation, therefore, we can 



compare theoretical results for LI = 200 nm and LI = 400 nm with the data in Fig. S7 1) and b), 
respectively. 



One can see in Fig. S7 that the nonlocal resistance is essentially independent of the injector 
length. The maximum signal is obtained when the detector is in the QSH insulator regime, as 
expected from previous nonlocal transport measurements in the QSH regime (see Ref. [9]). In 
the diffusive regime, we expect that a nonlocal resistance signal originating from spin effects 
would depend strongly on the system size. By a solution of the Poisson equation in the metal- 
lic T-structure we estimate that the stray diffusive signal for an Ohmic conductor doubles when 
the metallic leg of the injector is reduced by half. Therefore, a very weak dependence of the ex- 
perimental nonlocal resistance on the injector length excludes the possibility that our result is a 
diffusive stray signal and constitutes evidence that our samples are in the quasi-ballistic regime. 

We now provide further theoretical evidence that our signal originates solely from the presence 
of strong spin-orbit interaction effects, i.e. the Berry phase term (Sec. |XI]). Fig. |S8] shows the 



theoretically predicted nonlocal signal in the geometry of Fig. S6 As mentioned previously, this 



configuration is equivalent to that of Fig. |S7]by Onsager-Casimir reciprocity. First, the nonlocal 



resistance exhibits almost no size dependence as observed experimentally (Fig. S7). Since we 
have chosen the bulk impurity scattering time r to be larger than the boundary scattering time tqs, 
the agreement between the theoretical and experimental results constitutes strong evidence that 
our samples are in the quasi-ballistic regime. Second, the nonlocal resistance signal is negligible 
in the absence of the Berry phase effect. Since the Berry phase term is a direct consequence of the 



intrinsic spin-orbit interaction (Sec. XI), we conclude that the observed signal is chiefly due to the 
intrinsic SHE. However, as mentioned previously, we expect the Rashba spin-orbit interaction to 
also play an important role^"^, which could account for the discrepancy between our simple theory 
and the experiment in the magnitude of the signal. We believe that our supplementary experimental 
and theoretical results provide strong evidence that we have detected the helical nature of the QSH 
edge states via the SHE"^, and used the QSH helical edge states to convert the SHE into a charge 
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FIG. 1: Schematic layout of the two experiments on split-gated H-bar devices discussed in this paper. The 
dimensions indicated are those of the actual devices used in the experiments, a) shows the configuration 
where the current is injected into a metallic region (green, contacts 3 and 4). The spin-orbit interaction leads, 
through the spin Hall effect (SHE), to spin accumulation at the edges of the leg, as indicated schematically. 
The upper part of the structure (yellow) is in the QSH regime; the difference in chemical potential between 
the two spin states in the interface region is transferred by the helical edge channels to voltage contacts 1 and 
2. In b) the injector and detector regions are interchanged: the current is injected into the (spin-polarized) 
helical edge channels of the upper leg, causing a spin accumulation in the lower metallic region. The inverse 
spin Hall effect (SHE~^) converts the spin accumulation into a voltage signal. 
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FIG. 2: (a) An electron micrograph of the actual device structure (rotated by 90 degrees compared to Fig. 
1). (b) and (c) Gate voltage dependence of the longitudinal resistance Rxx (black) and Hall resistance Rxy 
(red) at B = 1 T, and the inferred carrier density, n, of a macroscopic Hall bar, 600 fim x 200 /im in size, 
fabricated from the same HgTe wafer as the nanostructures used in the experiments of Figs. 3 and 4. 
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Detector scan 




FIG. 3: Experimental nonlocal resistance data corresponding to the measurement configuration of Fig. 1 
a). In the bottom (green) panel, the gate on the current injection leg is swept, varying the area from p- to 
n-metallic conductance, while the detector (top) leg is kept in the middle of the QSH insulator regime. The 
red, blue and green arrows denote gate voltages where the injector region is p-type metallic, QSH insulating 
and n-type metallic, respectively. In the top panel, the gate in the detector area is varied at exactly these 
injector settings. 
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FIG. 4: Experimental nonlocal resistance data corresponding to the measurement configuration of Fig. 1 b). 
In the bottom (yellow) panel, the gate on the detection leg is swept, varying the area from p- to n-metallic 
conductance, while the injector (bottom) leg is kept in the middle of the QSH insulator regime. The red, 
blue and green arrows denote gate voltages where the detector region is p-type metallic, QSH insulating and 
n-type metallic, respectively. In the top panel, the gate in the injector area is varied at exactly these detector 
settings. 
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FIG. 5: Semiclassical Monte Carlo simulation of the nonlocal resistance signal in the setup of Fig. lb, as a 
function of carrier concentration in the metallic detector. 
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FIG. S6: Transmission coefficients for Landauer-Biittiker calculation of nonlocal resistance -Ri2,34, in the 
geometry corresponding to the QSH state as spin current injector, with the metallic state exhibiting the 
SHE~^ as a spin current detector (similar to Fig. lb of main text). 
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FIG. S7: Nonlocal resistance signal for injector in the p-type regime, as a function of gate voltage in 
the detector and injector length: a) short p-type injector (200 nm); b) long p-type injector (400 nm). By 
Onsager-Casimir reciprocity, for the detector in the QSH insulating regime iV*^,, ~ 0) this is equivalent to 



the configuration of Fig. S6 with LI = 200 nm in a) and LI = 400 nm in b). 
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FIG. S8: Effect of the detector size LI (geometry of Fig. [S6J equivalent by Onsager-Casimir reciprocity to 
Fig. [S7]) and Berry phase term on the nonlocal resistance. The dotted gray line denotes zero. The signal is 



essentially independent of LI, in agreement with experiment (Fig. S7 1. Furthermore, the signal is essentially 
zero in absence of the Berry phase term. 



